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REQUIREMENTS FOR IMPIjEMENTATION»OF KUSSNfflR AND WAGNER* 
INDICIAL UFT GROWTH FUNCTIONS INTO THE FLEXSTAB 
COMPUTER PROGRAM SYSTEM FOR USE- IN DYNAMIC LOADS ANA-LYSES 

il * « 4 « by RonEiild D. Miller asdjohn T. Rogers ' « • 4 * 

I 4 « k ^Boeing Coinioereial Airplane Company • « • • 4 * 


» • 4 


% • 4 


• 44 




»44 4«4 % 4 ^’ » 44 * 

. . » , , 1.0 SUMMARY* , . * g » 


* • 4 


k t 4 


k % > 


This document ^escribes the requirements«for implemeq^Uon of Kiissner and yfagfiei 
indicial lift growth functions iu^the NASA- Ames FLEXSTAB Computer -Program 
System (CPS) to represent unsteady aerod 3 rnamics, for use iji ^y^pmic gust, loads 
analyse:.. This stu^Jy was perfonjie: under N.^^ contract ^AS2-7729^ "Qe- jelopment qjf 
a FLEXSl^RComputer,jPrpgram” (ta^k|V, item I)^ ^ ^ ^ 4 4 « 

*♦ 44 ^4 4'i 4t k'4 

Use of the Kiissner apd^agner uimteady aerodynapaiq representations is rovi^wo^, the 
requirements for dynamic loads analyses are outlined, and the applicability to these 
requirements of the various program segments existing in the NA^A-Ames^ l.Ocfxx 
controls-flxed version of the FLEXSTAB CPS is determined. From this information, the^ 
modifications to the existing FLEXSTAB CPS required for creating a dynamic gust’ 

loads analysis capability are identified and discussed. ' 

- ■ 4 • » - « 4 ‘ 4 ' 

The conclusion is that a large nhmber of modifications and additions to the existing* 
NASA FLEXSTAB CPS woiild'he require before implementation of the Kflssnei* arid 
Wagner indicial* lift growth functions into* the FLEXSTAB CPS would provide* 
meaningful gust loads analysis capability.' 4 4* « 4 . « 4 • 


^ ^ 4 ' • ** » - 

^ ^ * 2.0 INTRODUCTION. . , ^ ^ 

. ' • *.* ■ ‘ •’ . 

* The NAS^ has developed the FLEXSTAB CPS (re^. *1) primarily^ for stabili^ ^aild 4 
[control Analyses of controls fixed'elastic flight vehicles using steady state ^eriodyaamvs ^ 
•and a ulo^-frequency^ approximation to unsteady aerod}mainics. Other qap^ilitie» 

^ include tke a|pili^ to caiculgt^ staltic load ^l^j^lAitioBS (LneVfia'^niLarr lo&ds)*an(^ti^:y • 
r contEO?r^<i?lire*iei 43 t(a: (lexiWe qirS-^; •#* 4 * * * 

t In order to4exp4na,tHe syst^^^ include dyu^ic loads. an^yses in volvmg %ctLve ^ 
control systenfs ai*d to^ma^e'it a4traetive for a variety of uses rangiifg from generalized* 

4 studies^ to actual desi^,4it istne^sriary to have ^he optionpf iming aerodynaimct theories ^ 
L ranging fr^m steady ^tSte to unsteady. Reference 2 exaiq^^d^ tfie feasjjbffity' of 4 

* : u„«.u 1:0. 
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I This si idy wa - liniiUd to examining only the use nf Kiissner and Wagner indiciai liPi 
^ growth unst< dy aerodynandc representati ns Spi-cific tasks were t( determine; 

j luirements for implementation and use of the KUssner and n'fif, ler indiciai 

I .Ti , . t.h jnti adv airodvnomit n, re; ntaiiuos »n tee iN'ASA-Ames FLEXSTAB 

; I.Uii.x tor dy uaruic loads anaivsvs 

' • ivdi itioiial program requirements to complete the preceding task 

t 

^ • Limitations and deficiencies which presently exist in the FLLXS ' Vt': i.OO.xx 
system for dynamic loads analyses 
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_3.0 GENERAL REQUIREMENTsToRJ)i^5/5^gi)jLDS 

ANALYSES, 

■» 'W T jf - •- «r -J - — - - • — - — 

To evaluate the FLEXS'^A^ CPS iu calciilatiijg dynamic gust loadai ii is first 

necessary to determine the general )^equi^ei9ent8,.of,d3{namu: loads analyse^ These 
requirements were established i^ r^ereye^ ^(^aje j'esjtatod^ye/or sub^^q^ei^ use ^ 


rwTv^ 


s « 


s s « « 


in this study. 

A general dynamic 'aeroe^astic^oads^a^alysm system ii^ludes tlm ^lloww^e^n^n^;, 


« # 

V 


1 >'Sr. #-#VS 

Definition of structvfral geometry and structural and aerodynamic panela — 

^ 4 » % flt u . ^ j- <r ^ ^ 

9. A'* « X i# • 
t?oi? (J 

* generalized mass and stiffness . ^ • 


4- ♦ * 


tl 


4. Definition of the aerodynamic model 


5. 


6. 


1 . 

• 

2. • Definition of the strficfur5l model^ 

3. * Calculation of structural vibration chafracteri^ics and determinatmi 

4 # « l» 

^ f ^ ^ f % • • • f I 

Definition of the control system model v ' 

Defyiition of excitation functions #« 

7. ^ Formulation of the equations of motion and load equations^" #4 % ♦ ^ 

8. Solution of these equations to determine: , *- 

^ ' t « tr ■ 4 * • 4 

a. Coordinate and load responses in either the^time.os,fnpqueqK:y<d«inain 9 asidt 

power spectral densitj^ (PSD) load parameters ; ^ ^ u 9 4 

4 ■ * -t ^ ^ •' 

b. * Roots of the characteristic equations 

« ^ . .5 - * i- - a «. * 1 # * » • < 

I The requirements and levels of sophistication for ea<A of these elements ^^eed to be , 

specified in order to estalblish a base for determining such limitations and daffcieQcies, 


4 4 « • * • f 4 

<4 m • 


' presently existing' in the FLESiSTAB CPS, as concern' dynamic loads ^aiyses 
capabilities. These general requirements are de^ribed in the following parag^ohs. 

» i • 

^ 3.1 GEOMETRY AND PANELING ^ * • 


* 




Sj^S. 

m. 


• • • 


Geometric >data describing the aircraft Components sfre re^ufreS before dynarnic 
aeroelastic loads analyses can be performed.* Panel chaYacterisHc^ pfer6nrat' to fhej 
' structural and aerodynamiiy model* representation must be* d'efi?ied; since* t''&e 
aerodynamic and structurait panels are seldom identical, separate ^paneling Schemes '* 
must be employed for ^ch. * 
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^ The structure may be ^(^es^nted by beam segment^ oi^ fl^it^ elemeijts for c^lcularimi * 
of either flexibility ^or s^ffness matrices If l^am sepnents are ixjodeled, ^the mafs 4 
representation usually copsists* of lumpe4 masses and Snertias-whi A, -whenf summed, ^ 
equal the total airplane c.g., weight, and rotary inertias. For a finite element structural J 
' representation, the mass is usually represented by lumped masses at the finite element 1 
nodes corresponding to the weight c.g.’s of small panels, and again total airplane c.g,, j 
weight, and rotary inertias are matched. The number of degrees of freedom resulting J 
from these structural and mass representations would be too lar^ for an efficient * 
dynamics solution. Consequently, the mass and stiffness or flexibility matrices are used^ 
to formulate an eigenvalue problem whose solution produces structural mode shapes^ 
and frequencies which are used to define a smaller number of generalized coordinates ^ 
than physical coordinates. Thus, the problem can be reduced to a workable size for 
dynamic loads analyses with little loss of accuracy (if enough modes are used). « ’ 

« . , % * ^ ^ 

, . , 3.3 AERODYNAMIC REPRESENTATION % i H 

» « . , ’ * * ^ ^ 
The aerodynamic* model is required to calculate the response and ejftH^tion^air forces | 

that are used to determine generalized air forces for the equations of hogtion 

development and to form air load contributions for the load equations. The degree of* 

sophistication required of the aerodynamic representation is dependent upon the tj^e pf^ 

dynamic aeroelastic loads analysis to be performed. For preliminary design and'gust^ 

loads, quasi-steady aerodynamics modified with Kiissner and Wagner indicial lift^ 

growth functions are generally satisfactory because of the large attenuation of the^gust^ 

forcing function at high frequencies. If mode stabilization, flutter suppression, ^r g 

stability augmentation systems are considered, then responses at intermediate or higher * 

frequencies are quite important, and a more exact unsteady aerodyneunic representation ^ 

must be used to correctly obtain phase and magnitude responses of the modes. % 

• ^ 1 

« , 3.4 CONTROL SYSTEM REPRESENTATION * ^ 

A control system definition is necesseiry tq assess the ef^cts «f aoiive %o^trhlSg 
responding to arbitrary inputs or feedback signals; e.g?, sitaBiliiy ’4u^mentat?oi% systqpa i 
(SAS) signals. The representation is dependent upon the t3rpe of d^amic aefoefei^ic^ 
analysis to be performed. In most cases, the control system cqn oe^'eprese^ted by linear^ 
systems amenable to classical control systeiy a^filysis ^d sjpt^e^is 'techniques' (rqf.^).|| 
In some fases (e.g., certain flutter suppressiomsystem^), no^liAeqf M^yms ^echi^q^s^ 
should be applied because of, elec^ical> comj^hent nonlinearitiesjl ser\% meehaniam 
saturation, or control surface movement limits. 

% % 

• * •3.5EXCltAp2>I^FijNCflONDBFINiTft)¥ ^ % % % 

• » • . * V • . • ^ • * • ' % ^ ^ 

Excitation functions may be* of saverat types. The^ Aa^ Anjisf ql" Sscjllito^jy 
arbitrary abrupt control surface motion er atmospheric turbulence. For dbntro> inputs, 
the control surface time history qr feedback control sig^ak serves to*, desceibe the' 
excitation function. Atmospheric* turbelence may be described liWtb dither*a discrete 
-nt di jendent or a continuous frequency-dependent moitel. In_^the forme^ a \raveform 


4 
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is sp^cifiedf; in tfie latitef^ a'pdWel' ^eclrum based on a statistical^e^^pt^oiT o^r^dom 
turbulence is defined. Acceptable models of atmoepbei^jc^tufbj^l^c^ ^re^ readily 


available (refs. 4 and 5). 




3.6 EQUATIONS OF MOTIONiAND LOAD EQUA^JON FORMTJLATTfON * * ' 

" • * « I 

Formulation of matrix coefficients ..for the equations of motion and load equations * 
requires structural, aerodynamic, and control system data as previously discussed. ^ 
Provision must be allowed for generation of either constant or frequency-dependent < 
matrices resulting from steady .state and unsteady aerodynamics, respective y, and 
either panel loads or loads at a reference location. There must be a 'capability “to* 
increase matrix size and insert, delete, or change individual' matrix’ eleraents»to> 
incorporate experimental data, additional degrees of freedom, etc. Load'equation raatrisf 
coefficient generation requires that load stations and an arbitrary load reference SRier 
I system be specified if shears, moments, and torsions are to be calcul£tted. In this casey , 
the appropriate panel inertia and aerodynamic forces must' be summed -tosthe k>a4 ( 
station and rotated into the load reference axis system. * 

^ 3.7 SOLUTION ROUTINES ••••■«•«««< 

• 4 - ♦ • - . 4 

Solution routines are necessary to solve the equations of motion and load equations Jm 4 

both the time domain (time history solutions) and frequency domain (in o^der to^ 
determine PSD load parameter and to perform time history solutions using Fourier 
transform method if unsteady aerodynamics are used). In addition.^ i^j i^ n^ces sary^to 
root the characteristic equation to determine stability characteristics. ^ ^ ^ j • * 1 


w ' r « 9 

The ability to perform dynamic loads analyses is dependent upon satisfying each of | 
these general requirements. The ability to .perform satisfactory ^srnpinic^oads analyses 
is determined by the level of technical sophistication gf each general requirement of the ‘ 
preceding elements. 

injLi 
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4.0 review’ OF THE INDfciALLlW GROWTH FUNCTION ' ‘ 
UNSTEADY SUBSONIC AERODYNAMIC REPRESENTATION 
^ FOR USE IN DYNAMIC LOADS ANALYSES 

The problem of representing unsteady aerodynamic flow assumed importance in Jhe 
‘ l930’s; subsequently, a number of aerodynamicists began presenting (solut^oi^ ^or 
various formulations of the problem. • * • * . 

Theodorsen first published in the United States a complete solution for the aerodynamic 
forces on a thin airfoil performing simple harmonic oscillations in a uniform 
two-dimensional incompressible flow (ref. 6). A solution for the aerodynamic forces on 
an airfoil subject to a step change iij angle of attack was developed by Wagner; the 
I result is reproduced in references 7 and 8. In these two cases, the solutiqn js expressed 
I as a product of the quasi-steady lift and a function known after the inve^igator’s name, 

^ which itself may be expressed in terms of Bessel functions. A transform relationship 
I exists between Theodorsen’s reduced frequeqpy-dependent function ^nd Wagner’s . 
i nondimensionalized time-dtependent function. * • I 

I The gust encounter problem has been approached similarly. Sears published a solution 
for the aerodynamic forces on a thin airfoil traversed by a sinusoidal gust (ref. 9j. A , 
solution for Kussner’s problem, which considers an airfoil encountering a sharp-edged 
gust, is reproduced in references 7 and 8. These solutions involve terms analogous to the j 


. Theodorsen and Wagner functions; the Sears and KSssner functions. 


The Wagner and Kilssner indicial lift growth functions are customarily designated # ^r) ^ 
and (Ti, respectively. Although these functions have a relatively simple form, they are ’ 
not expressible in terms of simple, well-known functions. Therefore, approximations ' 
written in simple algebraic terms were developed to facilitate use of the indicial * 
functions. The effect of finite span was considered- by Jones (ref. 10),. and additional ’ 
approximations of indicial functions were developed /or various finite aspec^ ratios * 
ranging from 3 ^ infinity. Compr^sibility effects ^were considered* hg ^a^l|ky and^ 
Drischle? (fefa. 11 through l3), ,ana %ov^e apprsKii^iate indicial «fugc|io|^ ^I^essimis^ 
were developed for various aubso^^c Milclvnumbers fjom 0.0 to (X7. • j m * * " ' 

'j • ' ''•ft 

A general procedure in past dynamic loads analyses for flight conditions iq the subsonic 
regime has been to use two-dimensional incompressible aerodynamic theory modified* to 
approximate finite span effects* (ref. f4); compressibility effects have been represented 
implicitly through \ise of compressibility factors (^)'to niodify,tho lif^cur^ slopes 
j^alculated using theoretical expressions for incompressible 4Iow. repres'enfelSon* of . 

^unsteady effects, tlie ^se of indicial function||has*been«x(en(le^ to include drbitrary> 
ai^oil motion, and ^gqst ' excitation functions through application of the superposition 
lintegraL This approach has been used extensively throughout the aircraift industry to 
• determine dynamic gust design loads for a number of sgbspnjc, relatively high as{^ct 
ja-atio, jet aircraft. ^Flfght test results have been*sbovim /o^^ee reasonably mel^ wit^ 
ftheor^tical predictions.** " 

* Richmond, L, D., "A Rational Method of Obtaining Three-Dimension^ Unsteady Aerodynamic 
Derivatives of Intersecting Airfoils in Subsonic Flow,” Boeing document D6-7401, 1962. * 

,^** Gilley, T. A. and Cast, R. D., "Theoretical and Experimental Frequency Response Functions for . 
the'B-52H Airplanes (WFT 1286),” Contract No. AF 34(601)-222^, 1966. 
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^ - f ; V m x 9 • - 

An analogous^ jft ^ 

^n^l^fn^ng a d^majpiif; oapobSt* iii tf 

«.- Jk • _® J\. i.^A.1 r 1 _ j. ^ j — ^ -m m 


«• 


( 

t 


>n|^ tj|fe Kef^iveMe Aa^or 

. ■ • • ^ 4 

aefodynamic IpnpuiptediBsftig^tMb fteVdy s^te VUC^i^i^||r^4e9iA 

AhA PLE^^TAB (^nj^inec^ vAthi tke IndtaStsKieous incide^e^^ ^n||icilIAifF 

growth ^£ctfons to ^velog aq^ approximate unstea^ Aerodynamic represent^ t%ma m P 

Unsteady Aerodynamic effects a^e approximated by considering lifti groijtl^t^oflpiar im 

A « » • • j 
. « A • » 4 


. accordance with indicial functions ^ ^e following forms: 

s 9 * * * * ’' ^ 

• Response indicial function approximation ^ 

» 0 m 

, , p * 1 - aj exp aj t - S2 exp «2 t 


I 




« I 


^ 2. ^Excitation indicia! function approximation 

• * 

• a « • • 

^ 'fee coefficients ap bi,|a|, and di are functions ^pect ratio and«sifosonic Mach’ 

^ ^ 1 Tn 1 *_ il „ 11 111 rm * 


* 4 1 


■ b j exp di t - b2 exp d2 t - b3 exp d3 t 


number. For supersoilic flow, these coefficients are generally assumed to be zero. 'The 
user is free to choose coefficients that tailor the indicial functions for compatibility with^ 


individual analysis requirements^ In theory, indicial functions can be simply developed 
~ from harmonic solutions; in practice^ this procedure requires ponsiderpble effort. • • * 

•'■■■ ' * A m 

• The actuations o^mi^tipn #is formulated irf reference *15 and refom^ul^tqi^ ijH vfareftc#2,4 
44 ^sing^tlfe ^asis o£ the approximate unsteady aei%dshi^ic represeiitaUoik cfosflrilfod' 


^sing’tlfe Tiaiis o^tlje approximate unsteady aei%d;j^m 

* previou^y, are: ' 

• • • • t • • ' 


rem^e|pq|a^o% desflrilfod* 

A 9 • • • * 
# * • • • 


• • (Strfcttr^ + SAS) + (Respond Aerodynamics) =(Gust Excitation Aerqdy^amiqs) • 

• • • ^ * • 

‘ [m H {q } + [m^ { q) + [ms] {q} + [m{| {q } ♦‘b + [ms] {q I *4> = |c3 jag*\p 


0 

0 «wbere • 


• * « 


• • 
• » 


M 2 ,^S ate fee appropriate structural matrix coefficients. 

# 


4 

• ■ * 

; • * • 
• 9 ” • 


« « 




, , « .. 0 m f 

M4, M5, and €3 are the appropriate aer^dynami® matrix coefficients formed u|^n|r 

^teady state aerodytiamics. ^ r ■ * * * 


.0 V 
» 


• • a 

q’s are generalized coordinates. 

1 « • • ' 

otg is tl)e gust angle (angle 
^fr^es^r^m, vqjooity with gust Histurbances) 

0 * 9 * * • - ■ 

0 ia Vfagper indicial lift 


♦ is t'he^Kus^er indicial lift 

* indicates convolution. 



a a A 


- - • 
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equations’ follow tiie same format as^^^atio'ns of motion. • ’ 

j Load= [Ml]|q} + [M2]{q} + [M3](q) + [M4]{q)*<l> +[M5]{q)*<I> + {C3} ag*^ 

1 where • • « . 

I # , * 

- - - . . , . , ■» * . , 

• ♦ Ml, M2, and M3 are load ^atrix coefficients at tha generalized coordinate 

f * • displacement, rate, and acceleration, respectively.-^ 

! ■ - - ■ ■ 

I * • M4 and M5 are load matrix coefficients of the generalized coordinate rate and 
f • • acceleration convoluted with the Wagner function _ 

C3 is^the load matrix coefficient of the excitation function convoluted with "the 
I * • Kussner^u^ction. * ®**-**tr^*,, ^ “ 

••• ^ 

The equations of motion, including- indicial- functions, may be integrated into the fdm* 
I (first-order differential equation) required in the existing FEEXSTAB CPS • by 
transforming the response indicial function into state form (see app. A). Although the 
matrix size increases significantly, this form allows a general application of the indicial 
functions. At one extreme, the same form of the indicial function may be applied to alj. 
, degrees of freedom or, at the opposite limit, each degree of freedom may be modified by, 
I a separate approximation of the indicial function as defined in appendix A. > * * 
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I*. :t ?<k : : : 

!•«•*••••••*••#*« 1 • • *«••«••*•«•• 

pi^ v^ri^u^ p(o^a|ps/)f,th| i)ff|fe§c^2^« 

^et|r^iii^e <^q^n^C|i Ic^d^ 

' ^er^djpqjnij r^pfesfntetion^ Iropi ^o^fr^^enfyvleQs^d^ti 

^ unsteady ^aerodynamics. The eval^troi^ was ^aged oq tlje g^e^r^l ^eq^ir^nqents ^et^iljd^ 

I in sectio^ to^e^ej with t^e Je^el of^techijical sop}iisticatipn^eg|u^e^ f(jr a (^i^nglc^ 
^ analysis. The basis for t^is^ stu(^ was tlia^ eyalriajiog, modi^ei^ t^ igclud^ i 

I Kussner and W^ner*in^icia^ lift growth unsteady |e^d^^y: ^epre|pi^a^oj^ 
addition, 'the evaluation *has iSeen expanded in several areas to furmer darifv tne^ 
» FLEXgTAB*limit^ti6nl. •** %.*.*«**»i*4**¥1ii4 

• •€■•••••■•••••! 

f rtgiire* I'idinftrils 'thi PLfeX^'rtVEP ftX^x^ ^o^afn^ c^^d^eS Ar^d^sftai? ftaWs ^ 
^a^pficaftidh, ’toother' wfth*thl efei5eifts1is^ed*iri*settiSn 3. * ••••••• 

!••»•«« • * I - «•« • «* 

5 . 1 «G 10 METini»A*J» PAItBimGP«SDf ••• ••••»« 

It* V >«*ff 

p ThetGeometrjt Eefinitioin ^iragram KJD) ^adcalaftes (fetiiled^difKithms %f*tllfe * 
t a«rodyaattic geometay «nd pianelingfi required in ther dew*s*rea(m*pi<og1‘anfts,*fr9n#tljfe * 
b generahiniput geoMetria data» Thi» aerodynamictpaneMng strtiaana w adeqfaate ioituA ik > 
Wpqrfomiina dynamic loads •nal 3 rsea. a «••••!•««»: «••••••! 

■ ••••••a****** «■•••• ••«««• 

^ « • • • • i.2iiSTmjQprvRAt»ni:iiRCSi;hP'ATiOiitfSiCaESifl)f» * • • • ^ 

^t^i^l|^tspc(u]^l Jn^u|p(^ (^^Tig^e^t Brqgram^L||^Cid ppea clasq^cf^ l|p£UP d^qpr|p 
y.o^oyn^q^ j^^stgfuq^i^alyin^dql qf qp ^t^c^ipl^o^pgnqpt^ ^egjr^^])|e^ 

^P|gr<^igga|^l^ (^ilj^i^iVt ||rit^ ^q ggag^lpa^b^q^ 

y:ompj^^ ^ j^oij^t^t^s^ffjge^ strgi^t Jiii^ ge^jjp^tg 
Qi^^u^oji ^ ^p^^nj^cl^b^ |p^ij|t^ 

ii'Wm 


4 

me'smicturar'geDraetry.'affa the lumpetr masses are ^nsidered to lie on this strai^tT 
fin^eft^c^i^.'rhfe 5aAr?oPtAs^epr^e^tanoS j^eftuae/ proper mcmeRng of a oSyj 
1 ^hSse^lfuftufal^i^s ft ifpi^dfit %nH (ft fflsela^e fo rfiofaP iiftrfia* fSc^, fom 9^ 

> #i# #e(?t Infts^hiAefri? afta^s8s.*Afth9i^ fti9 sRh9tiftaf flexftifft3FiAt9x^^ftsy 

1 fceftlfti<ftrft)<^fto#fci<Jlidf cftgfteftrtfrftdBnft cBn^UfttftrS df tfle ftoffy fti9lll 9rHs%)T 
I EliPffte#oi8 flre^ndt fc(ftu(ft(J*b^afts# <9 ifteft finftt^iifts^* Iftjg9'eF2,*cdin^Ai^ tftel 
I Asfte* dkdatl* 9LEX8TEB>0ES^Iisllc Ixll-Anftiet AaA fti#e^n«h(!lbifl, fll(fttlk#s * 

I tkeft Jbflftiafccfc^ ••••• •••••■•••••••ESAEttaaj 

I '%e«na 8 d|r^raBe«tMiAi Sf Mflki| 8 ^BrilcA ftoaidfestfoArdtaiy ^IHA ft *fttf tftsllnSl I 
|dagipeap%fii|e(||n 9 l^t^||t in liiftlftnding ftseAo^fifteftift. Voftidliift fte 9 t( 9 sftnH 1 
f i^#i4%e,sv * geqiiii^ft^iftdAeScl^ftnfledKittsftiittoava jftrft ' 

|a|d4Pft4t^ in 4PiP^^il^||j3|s ^enpendiculamtftijK i^attis wift iseao^ng 
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Figure 1.—FLEXSTAB Programs Reviewed for Dynamic Loads Analysis 
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Where: 

m = mass 

|0 = pitch moment of inertia about local elastic axis 
= roll moment of inertia about local elastic axis 
I|^= yaw moment of inertia about local elastic axis 
i = ith node 


Figure 2.— Slender Body Mass Representation 
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tlfc T(s*a^ l^rilbna^^dfejpevt l 

r ^ fllustrate^ ii^ figwe^^HfcweMer^ljgcjfuse in tiJe pr<%ran^to ^oW>e | 


Total mass 


Total pitch inertia 


Total roll inertia 


Total yaw inertia 


C.G. location 


^ ^41 Ai«cra^ ^«»rhich» the & locatioBa of thei^Kriiig«masaes «oi«ckle veiQT totals 
% * ^ location of the«tQtaViaircaafU canter gf #o4 tbe4>Qdy«cootnibai^ion 4p #i4| 

j ^ ^ i^]^aais«spj|^14 (^^plesv t]fi^#i4|g4iBOsaaact|p|ptiRCV»cq|if|fimatipq§)j^ 

k • 9 VI «rr(tn^U4ti)^al 

« kl 4^rya^ |^ayiq|[ ip^l^w^g^n^s^s^a^fl 4at^e,b^3^n]|is^s^ i^.^^ ^n|( 

I • 9 § large J)^v aj^d^r J)^V jJ'ud^ (a3{aj|pp^:^Cyi^rT|ft-jyo|| 

k • • V grr^jn^i^ tg^ «•«•«••« 1 

• *A*n* u^ •••!• ^ ■••/ ■• I* • • ? ^ ' 

A gexibk air^ane dynamic a^I^is r^u^e^ ^ ^ 

* Requirements, inclusion of me mas? cnara<^rtfltics, massunaonm^ of mertia at^e 


tt 
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reqi^re|n^t|^(^b^ ior e2||iv>% a^e^g|&r04ffoAcitcttafilol§ (# winf iftnftiil^ ^ 


I The m^s moment^of inerljia ^g) o^ tl^e^tqreds^a^lKeai^nd tbe«^ernudg9x \«iH14)e#oflfchi 

, sam? o?der £ Zmiy-z- fo^lightw^^t structj^rqa i^qjpi 9 s^«c(sirMat«wingtips. ClAsIk'di 
I this *exSmpl?, ne^lecti^ th| ineiji^ typjs 2Ijij wHaproducs* • 

'signlflcS^^e^oi^'in thTe l^n^in^ mo^e^t becai^e^f |h^nfr^£f%c1^ ■••••# 

■**’^*' ••«•••§ 


I The ^aftuTaffoi? ot' displacement* ^d rgt^icyis. o^ Fyoi^s,.oQ li/tigigB sp'||ce'^ 
I aero(9^%mic pahell is* based oh the defom^igp ^e^^r^1^t{U(^uypl|p(^eipciptiK>n 
■ the Mastic Sxife. *Tl?e "tw*is^ and displacement at^t^is^n^de \ptjj| the’ 

^ perpefcdftular*diita(fec#fAn?tIfe aerdd^amR ^nel centroid to^the ela^JiCj^^^js greiiq^t 1 
' to callulate? panel flefledticfhsr TTiis* approaJh Veglects the incjenaents ^n b^cjing^ qjiv.r' 

^ torsios deflection due* to #ie^ elastic *a^s 'diStsmce between the intersgc^oij| ^e = 

• perpendicular ^om tfie pattidl cfcifi-oM «hwf tHfe Itnictur^ node location (Ayjjf , 

' aa'W* 

' The Eiptamal Stmctural fciflue»ica ObeflficSeift |frol|rlin^^lC) is designed^to*tr|p^orm 
' stfuctyr^ matripe* oalouiat€<d ai^ng' finite "eRnftn^ fecfintqu*es, exterrmi to the 
( FLEXgTAB dJS*into stmctural matri«e#oP‘tHfe iftrA ftq8irdd*by*S^&feS^T^e*lM(S^ . 
^ body represent^tijgn Jiap_the*sa«i& deficiencies as those In'* ISICf that'ist tlie stracRiral 
f node points ara aasumqd to lie along the mean aerodynamic cei?terlii^.^K?s ?aiP r^ult 1 
( in largg Ijjc^izfd pffects op tbe«^tructure but the gross effects are smSll?' C6nsec(dedtiy, j 
this d^iciency,h^ qptinipq^l ^effects ftn the stability <derivati»e <cafculations fti * 

i but si^ifictyitly affects tjje generatiom of» the generalized inertia^and stiffn'es^ft^cH ’ 
, used in the (^nanuc equations o| mDtjyprv 

I Tlje use of cantilever, mpdes is de^rable for structural parameter studied aad^Copitibse 
; design cases where only a few stpuptura^ co,mponen|s need to be vaiaed (aeca4 ef *e#7). 

, System^^difications* required to gccommqdate j;autilever modes would dnaolweAni^nr 
changeAn^the iSid and ESIC programs because ^f t|Jieji^ure ^ tiiej-e#il<lang i»airi<*si 

In concTuAoh, SsfC*and/ to a lesser extent,^ E^Ip do . nqt jmeet* t||e apfj;u(|;urol 
represeiftaiSoii Aqiir^menfs as defined in s|ction3.. T^e^d^cjeqpi^ 41^1^ juntfbe'* 
j correcte# lftfdl'e*thes^ programs *are suitable for d^n^m^ j^oajjSgpq|il^^ 

I additiondor>mass*mbm'enf o? i8ei4^a*t«Tns ^r^fh slender |n(j^ t|jiiq,b(^i®s ^n|| (J) fo - 
slender Bodiesf tife fa^b/ftifjf H> ri&presen?both a curved elastic aiys not coincid^p^p^^^h *' 
’ the mea^anrodyBahiic* axis 'and lAasses Idca'terf off tAe elastic axis. T^eje j|'e(j|pijpV|^^|jUr 
' would alio ineBessitate*tHfe fcrfhufaftoif of flexfbiftty* matnces which ir^lv^e^d^^ioh^ 
and slop® caefficieihs tiufe (to momeht^. The capability is afso*n^ <^d eimer yiput 
\^ISf. rttodes • int<? ESlCr’’ dr fo’^rmulat^BTe^Plility matrices ii»fS}^Tr®fl> «h^ 

• cantileven aiodbsa camd>» calculated. Corrdbting' these deflbiehcies would req^uire a< 


^HyiUtiyc.# -it ''tf 9 9 0 ^ 4 


a •> A » 
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Where 


s 

i 


store 

wing struct mass 
0 



Figure 4.-Bending Moment Calculation by Force Summation Technique 




Figure 5.— Thin Body Modal Displacement Reference 
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^ ' The^ormal^modes prcjgranv (IjJVI) reqqjres »nput'fr*m tHfe s(tnictlral pro^am (I^Q. 

• MM^isesIlhis input J^o formicate an(i|Solv0 ag^ eigwalae |frobleSn to pr<^uce structural 

• * fttee-fi 4 e mode'^shapes and frequ^cies«thgit art> u#ed t# deflfie ■^eneftilized coordinates^ 
' ♦ and’ to calcJhate th^generaliz^ inertig and^stjjffnesit forces. Che stmfcturll iAput^ata, 

I . Snd* eigdh value problem fonnulaticjp and gpluti^n.^iare deficienl0 for <Eh<^ purj^osds of^ 
i ' HalcuRiting cahtilever^odes (ng^d body jpodea aj|p no^re^uiretf to (Saldlilat^c«(fttile^rl 

• Hiodes^ or inclndin^mdnien& of inertly in the jjiass maj^ix. Ponaequ 0 totly,*’redtrucl 5 iri 6 g 
4 and recodiag o^this program are necessary to sq^tisfyrtlj^ reqniremeiils ddfirrtd in^ 

•section l»for a satisfactory dyi^imic analysis. _ ^ >^# (0 ^ ^ 

• • # , - - ^ - - . # - * 

ff, ,0 » ► 5.4 AERODYNAMIC REPRESENTATION g\.IC,.SQ^SS)»^ 0 

*•« i*« ♦' * * _* ^ ^ « 4 0 * ' 

^The ^'unction of the ’aerodynamic influence coefficient grogr^ is tc^ calculate eiihe» 
0 subs9nic or supersonic steady-statfe aerodynamic influencp coeffiqjent (AIJl) mfitrices 
^ relating surface pressure to flow incidence, and to pro^ncq a iow-frequenoy 
^approximation to an unsteady AIC matrix relating surface pressure, to fldw incidence 
^ ^time rate of change.. The output consisting of elements of the ^IC jnaJ^ice^.is save^^for ( 
^ * Ujpe in the SD&S^ program where the-'complete AI^! rSatrix is assembled froq^ t^jese , 
* elements. m t * * ' m r ^ U 

^ '* * ^ ^ 0 . * *•* •• • ' * * * 0 0\ 

* "rtie a^odyymic,th9ory i^ FJjEXSTABwousidefs lineai? potent^l flow in^ot^ ^b^niC| 

* ^ an<^ supersonic ^gimps. ^Solu^ioqs i«e.« writteu as ■Integra# equSti^s invplving^t^e^ 

1 ^ ‘ strAgths of^ f^w sijagjjlarities^distr^uted orer aferodynaftic* mesifi sui?acfs ^nS pie|p 

* ' ifnes.^ Sfren^hs gf the flow ^singularity 'distributions ^re ifet^tmjned* using, 
I 9 afpproifimations based on^thosp used ip tJne fiadtq, etentont Aethod <fr Woodward. ^Thtug 
y # bo<^ and interference^ body n^ean pu^;faces*are represeiftedr by panels, ftich wiffe an 
^ Associated constant stren^^ vortex ^flow singularity.. Sfender bddytoieaff cetoter Aies are 
p * tKviddd into tine elements with flow singujarjties p-epreseHted by ^oub^fetrf whiJ^e 

f # strengths <fary qu^draticalTy.^ Thickn^s effec^ may .be represented by^softrce^ The 
a isolated ’brobleniS af e solved and th^ interf^encq in^idenae arising from these ^oliftion^ 
I calfiulat#d. TEis^intefference 'incidence ^s suppi^ssec^ w^h a . vortex distrihudion ftid*- 
r 0 isolated worteS ancTsohrce distributions. Filially, the pjessure is calenlated ahid tHI .^C 
is ^vmlabte from tMfe relStioh beween pressure and surface flow incidencer . ' # ^ 


I The steady state aerodynamic representation may be corrected with expefimantaWata | 
Fusing one of^seVerar schemes ^available in FLEXSTAB. In practice^ ^^a’nig^u^ | 
r corrections Uf t^e AIC- matrix are difficult to achieve, and the application af cerreato^ 
Lschemes ha^proven generally unsuccessful in the past (ref. 16 ). * m *" • *3 


The* low-frequ£ncy, approximation to uEBte 4 dj aerodygamjcs ij^ed^n ^LE^STAB hasj 
.severe firequencj* limitations suitable only for calcirtati(hi^f‘d^nhqp^| 

stabiFil^ty .^eVj^l't^es; FI^X^TAB oompPtes' fj-d'queh^-^i^epgnd^t dyn»mi*| 
^ coprespouding'dtQ Ihmt^ a%^-K),^etap T]|e iidienteit teestiictiAis 

■Areijfcen^ iSi d %n ^ agect^ratfo Ip Mach*nuiift)ei 4 ifeBitto%shi 5 

~ M^»fpr'suB 6 qplt,,ah^k<^l,^<<(M?.-*l)/M^ for 'sypersoniq. flow) , 
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"Result in an AIC matrix which is not suitable for dynarriic loads'* aif^ys?s 'invofvnig' ^ 
appreciable structural mode response. The envelope of these restrictions for a wing 
I having an aspect ratio of 8. 6- lB-521 in subsonic flow is shown in figurd 6.*Two flight . 

I conditions from -the B-52 CAMS s^u4y ^e^ superimposed on this figure to^slvpw their ^ 

I relationship to the testrictiona. *' »/, * \ m 

T ie use of steady state Alp’s and J,he instantan’eofis ‘angle of; incidence niodi^ed with 
t‘ e \Vagner indicial lift grpwth function, and thd* gust ’tingle modified with the^iissner 
i ii di lal lift growth function to represent uiisteadv aerodynamics has been the accepted ' 

; n etiiod used by the aircraft industry for performing dynamic gust loads amaly§es (refs." 

I 7 fc 17l. The aerjjdynamic influence coefficient program is .suiiablei fpr ^eneratinA 
steady state AIC’s for use in formulating the generalized respon^e^and gusff 
at rodynamic forces in a dynamic analysis. The low-frequency aerodynamii^ 
a pi ixiqjation is useful' for cafculhting the airplane static and dynamic stability® 

I a. ri\.atiye^ used in the dynamic equations of motion. It is important that the^^ ' 

I io>v-ifequency aerodyriamic approximation is not used. in formulating the generalized® 

I a> codynamic forces for ttfe elastic? modes (fig. 6), and especially for the gus*t generalized^; 


' in quency limits of fhe low-frequency aerodynamic. apprnxmation theory. #| 

i ‘ 5.5 CONTROL SYSTliM DijlFIN^TI^ON (SD&SS) ‘ , 

! Tl’.e control system model defined in the SD&rfS program can include aflerons. rudders, 
and elevators. However, if more than one control surface of'each type is defined, they 
must be interdependent; i.e., inboard and outboard ailerons cannot operate 
n'lependently of each other. Other control .<ar<'aces iduch as spoilers, tabs, or other 
! ,ir' iv< Controls such as side force generators) annot be’ modeled. Most present aircraft 
' >p- ra:e with some form of stability augmentation system tSAS); future aircraft will j 
I prf baoh have more sophisticated SAS systems functioning via a wide- variety of control J 
^ surfai es Any program capable of performing dynamic loads analyses must allow for a^ 
> variet V of active control systems. « 5 

* tv., • 1 

I 5.6 EIXCITATION FUNCTION DEFINITION (Sll&Sij.jm) ^ 1 

I « ’ t , , • ' J 

^The e.\ciiation function is defined in SD&SS as a time-dependent aerodynamic force at € 
§e«h- aei oci ntioid. These forces are ' useful ■ in« the equaUons pf potion formulatioi! 4 
^pqpc§^ur. to form generalized excitation forces. Time-dependent qhapta of the excitafio® « 
^angl§ of incidence defined" in TH include (*1 - cos), i«ti^,jai|d gquaie waves. To include i 
V, unsteady elTects, these instantaneous angles of incidqipce need to be modified with the q 
Kilssper indicial lift growth functions. » • * ■ 


5.7 EQUATIONS OF MOTION ..vND 
LOAD EQUATIONS bOKMLi.ATION (SD&SS) 


?ationp of motion are formulated ui the SD&SS program usir% data from the 
me defining the geometry, structural, aerodynamic and control .system 
n* ,n« There is no provision for including the SAS definitions in the 
ns S ice SD&SS is dependent on data from the n^evious nrogram.s, the degree of 
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Figure 6.— Reduced Frequency Envelope for the B-52 and for FLEXSTAB Subsonic Aerodynamics 
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piophisti^tiont o&tlie oPHiaF oT eacff ^ 

* 'programs. .PresenUy^ t^erp capabifity to formulate dyniamic load ei^ations;^ 

pwitheubthis ^ppb 4 lit^,^^F^EXSTA£ CPS cannot be used for dynamic loads Ctnalyaes.*^ 

UThe yis% a^trun^atgd jfiuml^r ^of modes plus the residual "flexibility of the remaining! 

Umod^, ^nclu4ii\g £^roel|istic effects as fonnufated ’in SD&SS* oP the FLEXSllAB J 
IlCPSjis not confidier^ a sati^faptory approach for use in dynamic load analysis. lUca«| 
^P|o<4ic^ lyis^ti^aptory resulte and its use is n’ot recommended. Th# reasons for -this ( 
L recomn^eiYlatioD sfe discussed in more detail in section 6. ^ f * 4 t ^ 4 

» • * « « * ‘ ' - t . 4 a 4. I 

. ^ , « . t 5.8 SOLUTION ROUTINES (CEIb, TH>f ^ ♦***:-* .* i 

^e* Characteristic Equation Rooting program *(CER) is a subprogram .of SD&.SSi;agdl 
solves for roots of the characteristic 'equations obtained from the equations of motion. 1 
The Time History program (TH) solves the equations of motion in the time domaip ^r a 
several forcing function shapes such as (1 - oos), step, or square waves. N 9 provision^i J 
i presently aifailable for solution' of “the equations in the frecuiency domain. . 3 

, . . , . : * : : 1 

f Since present ' design ’^st ^loads: criteria require the determinsj^op of Joads due to« 

F continVious turbulence, a solution routine must be qdded to solve the equations. o^ 
^ motion and' load equaticms using random harmonic analysis techniques. It must*be^ 
‘Capable of calculating steady state solutions for constant coefficient, linear second-ordei^ 
,^ifferentiah equations.^ From such steady s^te solutions, structural dynamic loacfe due* 
4lo sinusoidal forcing functions may be deterpained and statistical characteristics' ofload^l 
idue to conl^inuoos d^iucbulence calculated. The equations of motion must be solvbd*for f 
i tcoorcUnate i'reqqenqy response functions which are retained for use in calculating loa^ j 
ttranrfer* functions (simply load poefflcient matrices multiplied by coordinate responsea ( 
riat»eaeh^rqqi^qfy)^ L«ad power spectra can then be obtained from the product H>f "th# 4 
^Igust ^p 9 ct(a .,pnyd ^he.lopd tranrfef functions Squared! F*rom each load sp^trura, toot ) 
^mqaiv squq^e,lo^roo( nman |quare gust velocity (A) an(f number oPzero'crossinge 
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"residue^ wa^flfct' investigated ^ 

i60’s (reC 18). ^i#ca then, thi^ technique* for of 4 

^ugy 


The formulation Snd use of 
Sthwendler in the early 1960’s (ref. 18). ^ipca then, thi^ technique” for of 4 

sL.ibility derivatives and “for flutter calculations has been used successhiUy^I^v^viir,^ 
thi' use of thisrteehrtiqiie for flutter, contrc^l system analyses, and especially for 
loads anaJy§is,ican lead to sf^ificant errojsjf i%not exercised: * ^ 4 


i,«can lead to sf^ificant errors if care i%not ex 

^ m 9 • ^ 

1 o_i r ^ -.i»L ?j.' ^ ] 


*, m * * 


^e j>artitioned into two‘^efs:'*the ]o' 

f’ iuency explicitly retained modes Sn3 tlie higher frequency residual liiodes.' TheXa|^C| 

(J mo ■ 


Ir this techni^u*, 4he« free" vibration modei 


I 


a*- umption is t|h^ the inertia and dlimping forces for the higher frequency set of moc^s^ 
an much smaller t]pian the corresponding stiffnegs qndt thus* cafti he “n^lected; 

i.f, mq = cq« l^q^This can be considered as equivalent tq the genei'afiz^ ?na 88 ^ 
I b( oming zerq for those modes whose residual flexibility effects only are' ^e 
I considered (ref. 18). Figure 7 shows these effect^ ^n,j the% response' of *a^ 
I singla-degree-of-free^ojn lystem and the errow that can occuij, ^ 

1 ,«.*Bd**^''*J* mm'§rn4 

I For this example, residual effe^tanaw valM whhn^he forcing frequencies are«iac]Pl3k^ 
I than the actual resonance f;i;pqiie||C 3 i|o 6 the mode whose residufd Qej^bmty eifeet^dly^ 
I are included. Fo^ a multidegrqp-Qf-freedom system, where mod^l ^\^^,n§| OQCivs#tHlsl 
j region is rather fll-defined, Jbut the frequency ^f interest* must be n^c^ les^ than the | 

; resonance frecfuehcy of tiie last mode expKcitly* included. *MacNea| Sjchweadlarl 

1 recommended that the modes ex^ic|tly include l^hbuld ^e^ll' those in the^rqqi|enier4 
range from zero to the maximup frequency of interest ‘plus the next highei^ mode.f 

betweeg q^d^s is | 

* ~ • I 

• « ■ # 4 

n^lyses iBuat be «oPved in both the^tipe^a^jd ^ 
^ frequency domains. If ti<ne*hfetory*^s5u^ojjs ^re obtained, all the' rhodes emlmi^ 

^ included are esKiteds ^ the results qchiexedM}^ idfclidiffg^el^^a^^^^ 
* dependent upon tbet«n*uAt %iat each explicitlj^ ioclodad mo<le Ha^btfe^ exci*"^ — 



IS 

' the 


. — * 1 

degree of influence 'th0 dynamically*' deleted modes have on these retaihea mode^.^ 

fthe frequency of a -significantly' sfrong mode of vibration is close to the naturj 
^ frequencies of the modes whose residual effects onlj ye included, significant eir^s m 
be introduced because residual Stiffness-foiftes are included^withqutatbs rttetfltftg^h^ 
f forces to partially offset them. • 


r 


k • 


• > 


To show the effect^ of residual •flexibility* aSd truncated ,modes appCoSii^a^f^^S 
I dynamic analyses^solved ip the frequency domain, a two-degree-of-freedom mass-spnrife ® 
^problem without d^ny)ing was modeled* (fig. 8) and solved for the exact, tlfe ^ 
apprhxftiate residual flexibility, and tiie appibximate trunca 


I* In cpnjp n^^ ^ 
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Ki = K 2 = K 3 = K 
mj = m 2 = m 
Fj(t) = a j sin cot 




cjj^ = K/m 
co 2 ^ = 3K/m = 3coj^ 


1 1 

1 -1 


L whic^ coiwespoiHis 4b th< 

r # ♦ -» .♦• 

W 

I m m ^ ^ ^ 

' «u-iS 


Th» tw«P'phy»ic^®i|plffc^^pfln^*alieu?ated Djik5^-®#lj ujjire^* ^ 

iMk iikb ^^ir«i/«4uiC2/\1i'k4byMn • 4^ 4tR 4^ 


' m'^ * B»act)6dvrtlion:# ^ ^ 


I 



1 -(co/coj)'^ 3 -(co/coj)^ 


sin cot 


sin cot 




2K 


ai "1 


X, = 


‘1 

2K 


1 - (co/coj) 


2^7 


ai-2 


X. = 


1 -(co/coj) 


2 3 


sin cot 


sin cot 
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ai 1 

2K 


Xo = 




1 

1 -(co/ojj)^ 


1 -(cj/cjj)^ 


sin cjt 


sin cot 


rAlTiiou^ the ^eq^Suzed inertiff and damping forces are Mnsi^ere^mSl 
* assumptions and are neglected for the modes whose residual effe^'only are inclined, 

^ this does not mean that the accelerations of these modes are zeroi butf*only tl^^*they are^ 
I small. Thus, to, obt^h accelerations, the displacement solutioq§ can be differentiated ' 
[ twice, yielding thohfollowing: *'•' _ * 

r ^ ^ *' * . 

Xfp ■* ^ 


m « * 

i m Exact Solj^tioi^ " » ~ 


2K y _ 
^1 ’ 

2K 




o? 


l-(co/u)j)2 3-(cj/coj)2 


Xo = 


CO 


CO 


1- (co/coj)2 3-(co/coj)2 


sin cot 


sin cot 


Resjcfual Flexibility Solution: 


2K‘ 


^1 = 


2 2' 
(jJ^ co'^ 

^2 3 


1 -(co/coj)-^ 


sin cot 


2K 


Xo = 


CO 


1 -(co/coj)2 


CO" 

3 


sin CO t 
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^ Modal Truncation: (U 2 “ I 


2K 


ai ‘1 

2K 


X, = 


‘1 


Xo = 



sin cot 


sin cot 


r . 

t 


These solutions are sho^ graphically ip figures 9 and 10. When the model includes 
viscous damping ,*>the adution becomes considerably more complicated. This problem 
was formulated apd solved 01 ^ a 'digital comput^/, andsthe results for the displacement 
and acceleratien of mass I are sh^Wq^n figuoe ^ ^ , 3 , 

The results shown in figures 9 through 11 illustrate the fo^owing jaunty ^ ^ 

.♦ e ^ I 



the 
included 
solution 
resonance 


p iii calculating displacegi^ts. Ih c^ '^prdfvide (^nsid 2 >ral>le iiflorovemen^ &T M 

\ « truncated-moiJe'^lutiolB. - 

f r 

Tor frefluendes below t^xe fi*t ^sonance^freque|^(^, iaclqd^lf the^stal^cf ^d|>onse,^lT^^ ♦ « 

I residual«flexibility approximalioi» produce# a elisQl^|nienteai^*accelei»ti(m^^o4er^ ^ ^ 

exact* solution* than the t^ndbte^ iy> 8 fes *solutiop.^ H«ve\jpr,^vrtiei» d^p.^n 2 * i#» ^ . 
ided in the sy^;^, the res^jll^ 'flexibility solution be^^s^ ^p^t^frota tl^e ^act^ 

'in at a lower frequency thpn without damping. At ffeq^Aicips e^ve fhe «ecfnS^.^^ f 

nee fi^quency (co^ = V 'displacement of*botb exact <an<f ftuiftat%ii ^ * 

modal'* analy|is' approaches zero, wl^^ms' the residual flexibflity^ amxrqpc^fAioik 
^ ' approac&es% constant. Similarly,' the accmq):at|ons for the exact anc|^tl^neatad ' 

^ analyses ^^proarh .^onsthnta, ^ereats Ibe .acceleration of the aes^uiy ^ # 

approximation div^rgps a square of die frequenpy. Hbwever, the Itsi^ira^ lfbxi^ni||r * 

[ solutions are invalid atffre^uendles wh^ ^^oalth qr^e^ceed the,n^tuAil ^c^e^ilb ^ 
p • the modes not explicitly enci^d^.* ^ »*• ,'^♦11 ,.**.** 

•*> » ' ^ * ’ *# ***^«* 

^Fdr random harmonic analysis tp<3m^ueb, thbonetic^ derivatioi^ F^direa ^ * 1 

calcnilation of polutionsdorfc^cii% faeqi^hcies ranging froAi zero to^nfihit|K i^l^acflfc^ ^ ' 
>thq u^per limit of the freq^ti^*¥a 4 ^e ^ ^ome flni^ ^Ide. fllii^ ireq||e^]|Ks, ^ * 

^ .asu^lv^s^b^i^e d^b a^pd 09 the* joteniit yeiy ^)%_Jp«x^[i^^frfn(^^a|i4pii;iq|| 


4t 


'^^efng,.lIt&Q were used to fgri^l^te^aiid wive ^a- two-degree<of^eed% modeh with 

viscous danping. 1' •“ ^ - 

■ » . -W ^ 
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Figure 11.— Response fora Two-Degree-of-Freedom System With Viscous Damping 
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f requency“anf aTSP' convergence ' JKe ?Sl) paramet^ (for ^ample, 'Atmospheric 

turbulenc^ decreases as appB 0 ximate^*(l/freq) 2 ). (Jpnsequently, niode^^^^ose 
frequencies are slightly diigher ,JJian this finite •cutoff v^lue, ^s“*well as thos# whose 
frequencies" fall 'm this range, must be included to prevant significant errors froiy^beuig 
introduced into the anal 3 (sis ^en residual flexibility approximations are induced. A3 
can be seen from'figures 9, 10, and 11, the error introduced by using ‘Fesidual flexibili^ 
approximations is much larger at frequencies above fiie reasonance fre'quencyj^anjthat 


( error if the mode is left out, as in truncated modal, analyses. Table 1 suram^i^^he 
results’of the precedtng exan^le'as the forcing frequency approache’s infinitjii. i‘ ^ 

im ^ 


Table 1.— Response of an Undamped Two-Degree-of-Freedom System 
as the Forcing Function Frequency Approaches Infinity 


Solution 


Residual flexibility 


Modal truncation 


I is in the iqodeliqg of all elasti<^ffe<^s in thp rigid Ibody re^onse of the aircraft. 4 
i ltowever,''*since i^is necessary to have- modes explicitly included whose fr|gu%hcie8 . 

range from zero to a level higher than, the maximum frequency of interest which is' 

' itself fairly high (order of 15 ter 20 Hz), the increased complexity of iqplud^g’^residual 4 
f efTects,-Tather than using truncated modes, is not jus'tified. Because of th^ iiifcreadfed 
f complexity eind uncertainty at what frequencies the’ solutions become invalid, iWia * 
H-ecommended that this technique not be employed for dynamic loads analyses, and that 
^only truncated modeafTie used>with.a sufficient number of modes included to'encodipass , 
the highest frequency of interest in order to have reasopably accurate resuUiR In j 
’ addition, the modes considered must include* at least one.^ending and one^orsii^i Mode 
;foiv each lifting svffface. This is. required to obtain the plastic effects of*all 4^t^ices, ^ 
I including cases where the structure is d 3 mamically stiff but staticatly flexible. *** 
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F(t)/6K 

OO 

OO 

0 

0 

F{t)/2K 

F(t)2K 



_l 8-1/Z X 11 INCH CROPS 


I « 

U 

1 


. ' i 

The following modifications or ‘ additions ^o’ the |’LEXSTAB Gpll are^ n^c^sarjfc to 
provide the capability to adequa^l3|j)erfon^d]^amic Idads ^a^^sqs. ^ ^ ^ ^ 

ic axis r^^resentati^\|o allovi^a 4!u^ed, 

* \ ^ ^ V 4 ^ ^ 4 

♦ V' Ik'* 

ude mass ^c)|r<%ap' ine|1;ils 

.v' *'v»' 

_e prop^ mass disfcbution» aad ^ 
incl^d^ bending ^ ^ 


Modify the slender body elastic 
located elastic axis. 

% 

Modify the slender body mass fepresenta^ion to include 
off the elastic axis. ^ i| I ^ % 

> ^ ^ V ‘ . % ^ \ 

Modify the thin body mas^t representation to include 
rotary inertias. ^ 

, , t ' . » I V 

Modify the aerodynamic reference point interpolation rout 
^nd rotation slope changes between structural nodes. ^ 

^ I ‘ ■ ‘ V - « 4 4 , 

Modify the flexibility matrix formulation to allow escalation qf (fieptacen^h!^ and ^ | 
rotations due (o moments idiadditiem to current wplaqemei^? and r^a^om^u^t^l - , 
force. 4 ^ ‘ . V' ^ V' 

Incluck tke capability of using cantilever toode^ ii^tife an^yssil ^ 

\ ^ ^ ^ ^ ^ ^ ^ % ■ 

Modif^the program to calculate either cantilever ^ rree-fr»e4iOTes. |a ■ w ^ 4 * 

,• 4*' 

b^clhde the capability to calcula^ |oau equations coimis^n| of ne^^* 

^ momei^s,k torsion, net panel forces, l^ellKiynanMc^i^el for^s^ % 

^ vl^locities^i^ displaceiuents. k 

-4* ji 4.^1 4.^ * k “*1 ^ 

.Add thq d|)^illty^to inllime th% representrtio^ i|| 4^4<ma ^o^plFmoko^ 

1L^ A ^ ^ ^ 

Add a nev^oiution routine ta sdive the ^quarions of lyo^n andMoi^ a ^ 

the fSequencjkdomain, - .% li. ^ 
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Implementation of these momficatidns or Additions to tlv FLEXSTAB^CPS'will make it • 
suitable for dynamic gust loads prqvidecf that«,if a SAS system is included, it primarily 
I influences only the* low-frequency rigid body airplane response. However^ the Kiissner 
• and Wagner indicial lift growth. function representation of unsteady aerodynamics is not 
satisfactory for analysis of aircraft^here mode suppression or flutter suppression with 
! active controls is required. This type of analysis requires accurate phase and magnitude 
i relationships between modes, which requires a more exact unsteady aerodynamic and 
I structural representation as outlined in reference 2. The advantage of the analysis 
system outlined 4n reference 2 over this system is that it is satisfactory for dynamic 


gust loads analyses employing techniques such as the Kiissner and Wagner unste£idy 
aerodynamic approach and for both dynamic gust loads and active, control analysis 
using more exact theories. It thus -appears that implementation of fhe system proposed * 
in reference 2 is a much more effective |ipproach than Would jesult from implenMnj^tioi? 
of the required changes determihed f^m thiaistudy. * * 

Boeing Commercial Airplane Company * ^ ^ ^ ^ * 

P.O. Box3707 ^ •* 

Seattle, Washington 98124; Septembg^ 1975,^ * • ^ ^ 



30 


8-1/2 X 11 INCH CROPS 


' * • . APPENDIX A ♦ 

I ^ ^ * ** ^ ^ 

i EQUATIONS OF MOTIOIg INCLUDING THE WAOT^ER INDICIAL 
LIFT GROWTH FI 4 KCTION IN STATE ^ORM 

i _ A* 1 ^ * 


# 

’ ■* 


I When an unsteady aQ^odynamic representation losing indiciai lift |r6wtJPfunc^ons is 
, consider,ed, the dynapiic equHions qf motion fomin aircrWt may b^writteiiiin t^ form: 

r 1^. . li. . ? It - * » ' <• • _ 

[Mjjq} + [M 2 j{q} + [M 3 ]{q} + [M 4 ]|q|*(/. + [MjJjq} *0 = {C3}dg*0 (A-1) 



L5mc?3 


q(o) = q(o) = O(g(o) = 0 


£[q(t)] jC[F(t)] = jC[/^q(r)F(t-r)d7] 
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19 ^ 

* Then the Laplace transform of'equation (A-1) iszj. 


:|JI 


[Mi +sM 2 + s2m3 + (M 4 + sM 5 )s • 0(s)]{q(s)} = {C3}-ag(s) • s ' i//(s) (A-2) ^ 


!■ 


where 


etc. ' ^ 


4 ’ ~ ^ 






^ Consider an approj^imation to the Wagner indicial lift growth function having the form: f 

Jfc.- ^ ^ ^ 


k 


-ait -^it 
0 = aj - b jC - cje 


I then j 


I 0(s) = aj/s-bj/(s + aj)-ci/(s + ^l) 

i 

S- 0(s)= ai -s • bi/(s + ai)-s • ci/(s + 0l) 

^eflne vectors r to write the Wagner mdidarfunction'in state ?5rm 
_ ** ^ ^ aj{rj(s)} = bj/(s + aj){q(s)| 

ai{r 2 (s)}= cj/(s + /3i){q(s)} 


1 


(A-3) 


I Substitv 

I • 

» » I 


Substitutie equation (A-3) into equation (A-2); then equations describing the system are: 

[Mj t- sM 2 + s^M 3 j|q(s)|+ aj • |^M 4 + sM 5 j|q(s) -?j(s) -? 2 (s)| ^ (A-4)^ 


K «• 

1 r- 


•1 


= {Fi I* ag(s) • s • 0(s) 
s • bi|q(s)| = aj • (s + aj){rj(s)| 
s • cj |q(s) }= aj • (s + (3i)^2^®)} 


J 


[la-® 
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